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Summary 

Holographic  interferometry  and  computer- 
assisted  tomography  (CAT)  are  used  to  determine 
the  transonic  flow  field  of  a  model  rotor  blade  in 
hover.  A  pulsed  ruby  laser  recorded  40  Interfero- 
grams  with  a  61  cm-dlam  view  field  near  the  model 
rotor-blade  tip  operating  at  a  tip  Mach  number  of 
0.90.  After  digitizing  the  Interferograms  and 
extracting  fringe-order  functions,  the  data  are 
transferred  to  a  CAT  code.  The  CAT  code  then 
calculates  pressure  coefficients  In  several  planes 
above  the  blade  surface.  The  values  from  the 
holography-CAT  method  compare  favorably  with  pre¬ 
viously  obtained  numerical  computations  and  laser 
velocimeter  measurements  at  most  locations  near 
the  blade  tip.  The  results  demonstrate  the  tech¬ 
nique's  potential  for  three-dimensional  transonic 
rotor  flow  studies. 

1 .  INTRODUCTION 

On  many  helicopters,  the  rotor  blade's 
advancing  tip  encounters  transonic  flow  during 
forward  flight.  At  these  high  Mach  numbers,  the 
rotor  blade's  performance  suffers  from  compres¬ 
sibility  effects  that  often  cause  shock  waves  to 
form  near  the  blade  tip;  the  shocks  can  extend  to 
the  acoustic  far-fleld  (1).  Through  theoretical 
and  computational  investigations,  researchers 
attempt  to  understand  the  local  shock  generation 
of  hlgh-tlp-speed  rotors  and  Its  propagation  to 
the  far-fleld.  However,  because  of  the  problem's 
complexity  and  the  difficulty  of  obtaining 
detailed  experimental  information  about  the  flow, 
accurate  means  for  confirming  transonic  rotor- 
blade  designs  have  been  notably  lacking. 

In  an  attempt  to  describe  the  transonic  rotor 
flow  field  and  to  resolve  the  problems  associated 


with  it,  promising  theoretical  models  [2]  and 
numerical  codes  [3,4]  have  been  developed.  The 
numerical  codes  compare  favorably  with  blade- 
pressure  measurements  [5],  but  are  not  yet  veri¬ 
fied  at  points  away  from  the  blade's  surface. 

Holographic  interferometry  is  an  effective 
diagnostic  technique  for  making  transonic  flow 
measurements  [6].  Previous  Investigations  [7,8] 
in  which  two-dimensional  flows  over  airfoils  were 
studied  show  that  accurate  quantitative  informa¬ 
tion  is  obtainable  using  holographic  interferom¬ 
etry.  However,  the  transonic  flow  around  a  heli¬ 
copter  rotor  blade  is  three-dimensional  and 
requires  a  tomographic  technique  to  compute  the 
correct  flow  Information  from  several  Interfero¬ 
grams.  To  date,  most  applications  of  this  tech¬ 
nique  have  been  limited  to  axisynnetric  flow  or  to 
simple  three-dimensional  flow  with  a  small  model 
under  ideal  laboratory  conditions  [9,10]. 

This  paper  discusses  the  procedures  necessary 
to  obtain  quantitative  measurements  of  a  tran¬ 
sonic,  three-dimensional  flow  field  near  a  rotor- 
blade  tip,  using  holographic  interferometric  data 
and  computer-assisted  tomography  (CAT).  Though 
most  helicopter  rotor  problems  caused  by  shock 
waves  occur  during  forward  flight,  this  experiment 
Investigates  the  steady  problem  (hover),  which 
simulates  many  physical  phenomena  of  forward 
flight  [11].  The  method  for  recording  Interfero¬ 
grams  and  example  interferograms  is  included,  and 
the  steps  required  for  extracting  quantitative 
Information  from  the  interferograms  are  outlined. 
The  technique's  potential  for  measuring  three- 
dimensional  transonic  rotor  flows  Is  demonstrated, 
and  the  results  of  using  this  technique  are  com¬ 
pared  with  those  from  previously  performed 
numerical  computations  and  laser  velocimeter 
measurements. 


2.  BACKGROUND  CONCEPTS 

For  the  experiment  to  be  successful.  It  Is 
necessary  to  1)  record  high-quality  interferograms 
near  a  rotor-blade  tip  from  multiple  viewing 
angles,  and  2)  implement  a  suitable  CAT  code  with 
the  Interferogram  data.  Familiarity  with 
holography,  holographic  interferometry,  and 
computer-assisted  tomography  principles  provides 
the  necessary  insight  for  understanding  this 
technique. 

2. 1  Holography 

Holography  is  a  two-step  imaging  process  in 
which  diffracted  light  waves  are  recorded  and 
reconstructed  (12,131.  The  first  step  is  record¬ 
ing,  or  storing,  the  hologram.  This  is  accom¬ 
plished  by  dividing  a  single  coherent  laser  beam 
into  two  beams  and  exposing  a  photographic  film  to 
the  two  light  waves.  The  object  wave,  which  is 
the  wave  containing  the  flow  information,  passes 
through  the  measured  field  (the  air  near  the  blade 
tip  in  this  experiment).  The  second  wave,  the 
reference  wave,  passes  around  the  field.  By  add¬ 
ing  the  coherent  reference  wave  to  the  object 
wave,  the  photographic  film  records  a  high- 
frequency  interference  pattern.  Once  the  film  is 
developed,  it  is  known  as  a  hologram,  which  is  a 
complicated  diffraction  grating. 

The  second  step  in  holography  is  reconstruc¬ 
tion,  or  playing  back  the  hologram.  This  is 
accomplished  in  two  ways.  First,  a  reconstruction 
wave  identical  to  the  reference  wave  Illuminates 
the  hologram.  The  hologram  diffracts  the  recon¬ 
struction  wave  and  produces  a  replica  of  the  orig¬ 
inal  object  wave,  forming  the  original  object's 
virtual  image.  In  the  second  method  of  recon¬ 
struction,  a  reconstruction  wave  which  is  con¬ 
jugate  to  the  reference  wave,  illuminates  the 
hologram.  The  hologram  diffracts  the  conjugate 
I  wave  forming  the  original  object's  real  image. 

The  real  image  may  be  photographed  without  the  use 
of  a  lens  by  placing  a  sheet  of  photographic  film 
in  the  real  image  space. 

Several  Important  characteristics  of  holog¬ 
raphy  are  applicable  to  the  experiment  at  hand. 
There  are  very  few  geometrical  constraints  in  a 
I  holographic  optical  system;  thus,  holography  can 

be  applied  in  a  large-scale,  nonlaboratory 
environment.  Note  that  recording  and  recon¬ 
struction  of  the  hologram  can  be  done  in  different 
locations  if  the  reference  wave  is  reproducible. 
This  allows  the  reconstruction  to  be  done  in  a 
laboratory,  far  from  the  harsh  environment  in 
I  which  the  hologram  was  previously  recorded.  The 

reference  wave  serves  only  as  a  method  of  record¬ 
ing  and  reconstructing  the  object  wave.  Thus, 
a  hologram  does  not  produce  quantitative 
information  about  the  field  of  interest.  To 
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obtain  quantitative  information  (in  the  form  of 
interference  fringes)  an  interferogram  must  be 
recorded . 

2.2  Holographic  Interferometry 

Holographic  interferometry  is  the  interfero¬ 
metric  comparison  of  two  object  waves  recorded 
holographically  (see  ref.  12  for  further 
detail).  In  this  experiment,  the  two  object 
waves  are  recorded  sequentially  in  time  with 
double-exposure  holographic  interferometry.  The 
interferogram  is  recorded  by  first  exposing  a 
photographic  file  to  a  first  reference  wave  and  an 
"undisturbed"  object  wave.  Later  in  time,  the 
same  photographic  film  is  exposed  to  a  second 
reference  wave  and  to  a  second  "disturbed"  object 
wave. 

When  the  holographic  interferogram  is  recon¬ 
structed,  the  virtual  or  real  image  shows  the 
object  (the  transparent  field)  with  an  inter¬ 
ference  fringe  pattern.  The  fringe  pattern  repre¬ 
sents  the  difference  between  the  "undisturbed"  and 
"disturbed"  flow  states.  The  Irradiance  of  the 
reconstructed  wave  is  proportional  to 


where  and  Uq2  complex  amplitude  of 

the  first  and  second  object  waves,  which  can  also 
be  written 

(1)  I  s  2A^[I  ♦  cos(A*)] 

Equation  (1)  represents  the  interferogram  with  a 
fringe  pattern  of  dark  and  bright  fringes  of  con¬ 
stant  optical  path-length  difference  (OPD)  A^. 

AA  is  given  by 

(2)  AA  =J[n(x,y,z)  -  n^lds  s  NX 

where  n  is  the  refractive  index,  N  is  the 
fringe-order  number,  and  X  is  the  laser  wave¬ 
length.  To  determine  the  flow-field  properties, 
the  line  Integral  of  Eq.  (2)  must  be  inverted  and 
solved  for  n(x,y,z),  the  refractive  index  at  a 
specific  point  in  the  field. 

In  a  two-dimensional  flow  (i.e.,  the  flow 
over  a  fixed  airfoil  in  a  wind  tunnel),  the  evalu¬ 
ation  of  Eq.  (2)  is  simplified.  Since  the  refrac¬ 
tive  index  is  assumed  constant  along  the  probing 
rays,  equation  [2]  becomes 

(3)  n(x)  =  Oq  ♦  Nx/L 

in  a  horizontal  plane  (y  =  const)  at  any  location 
z  across  the  test  section  of  width  L.  In  a  two- 
dimensional  flow,  the  fringes  on  an  interferogram 
are  contours  of  constant  refractive  index,  thus 
the  refractive  index  at  any  point  in  the  field  can 
be  determined  from  a  single  interferogram.  In 


the  general  asynnetrlc  case  (i.e.,  the  flow  field 
near  a  rotor-blade  tip),  however,  equation  [31 
cannot  be  used.  A  set  of  Interferograms  must  be 
recorded  at  different  viewing  angles  around  the 
field  In  order  to  reconstruct  the  flow  field.  The 
refractive  Index  at  any  given  point  In  the  flow 
field  nay  then  be  obtained  by  applying  the  nethods 
of  conputer-assisted  tomography  to  Invert  the 
integral  equation  (2). 

2.3  Computer-Assisted  Tonography 

Tomography  is  a  mathematical  technique  for 
reconstructing  a  three-dimensional  field  from  Its 
two-dimensional  projections  (see  refs.  14  and  15 
for  a  wide  variety  of  applications).  A  projection 
of  a  three-dimensional  field  Is  the  fringe  pattern 
recorded  on  an  interferogram.  All  methods  require 
multidirectional  projection  data  of  the  field. 
Tomographic  codes  develop  in  two  directions: 

1)  Iterative  algebraic  reconstruction  techniques 
[16],  and  2)  Fourier  transform  techniques.  A 
version  of  the  latter  method,  termed  filtered 
back-projection  [17,18]  appears  most  suitable  for 
this  application. 

Most  Fourier  transform  techniques  employ 
back-projection.  Projection  data  from  the  field 
are  recorded  in  one  plane  at  several  aslmuthal 
angles  around  the  field.  For  example,  one  pro¬ 
jection  of  a  uniform  absorbing  disk  is  shown  In 
figure  la  (taken  from  ref.  19).  Beyond  the  disk 
boundary  (no  path  length  through  the  disk),  the 
light  ray’s  OPD  is  unchanged,  producing  no  inter¬ 
ference  fringes.  Near  the  disk  boundary  (short 
path  length  through  the  disk),  the  light  ray's  OPD 
is  changed  slightly,  producing  a  few  interference 
fringes.  And  near  the  disk  center  (long  path 
length  through  the  absorbing  disk),  the  light 
ray's  OPD  is  changed  substantially,  producing 
several  interference  fringes.  Similar  projections 
(fringe  number  vs.  position)  at  different  azi¬ 
muthal  angles  are  also  recorded.  Each  projection 
is  then  back-projected,  or  smeared  back  along  the 
direction  in  which  it  was  recorded  (fig.  lb). 
Values  are  added,  point  by  point,  to  form  a  recon¬ 
struction  of  the  field.  Unfortunately,  simple 
back-projection  produces  an  undesirable  spoke 
pattern  which  severely  degrades  the  quality  of  the 
reconstructed  field. 

To  eliminate  the  spoke  pattern,  the  back- 
projections  are  filtered.  A  one-dimensional  con¬ 
volution  (indicated  by  an  asterisk  in  the  follow¬ 
ing  equation)  is  performed  between  each  projection 
and  an  appropriate  filter  function  (see  ref.  20 
for  a  discussion  of  filter  functions)  before  back- 
projection,  as  shown  in  fig.  2a  (taken  from 
ref.  19).  Each  filtered  projection  is  then  back- 
projected  over  the  reconstruction  space  (see 
fig.  2b).  The  negative  side-lobes  Introduced  by 
the  filter  eliminate  the  spoke  pattern  during  the 


polnt-by-polnt  addition  process.  The  convoluted 
back-projection  is  given  as 

n(x,2)  =  fj-'  Ng(t)  •  F’’{|f|}«(x  sin  e 
♦  z  cos  e  -  5)d5  de 

where  n(x,z)  Is  the  reconstructed  field  function 
(refractive  index)  In  a  horizontal  plane 
(y  s  const),  and  Ng(5)  are  the  cross  sections 
through  the  fringe-order  function  at  angle  t  to 
be  convolved  with  the  inverse  Fourier  transform  of 
the  abs-functlon.  The  integration  of  the  con¬ 
volved  "projection"  is  done  along  lines  expressed 
as 

6(x  sin  e  ♦  z  cos  e  -  () 

and  overall  projections  0.  With  many  projec¬ 
tions,  this  technique  yields  an  accurate  recon¬ 
struction  of  the  original  field  in  one  plane. 

This  procedure  is  repeated  using  data  in  several 
planes  to  yield  a  reconstruction  of  the  entire 
three-dimensional  field. 

3.  pr(x:edure 

Several  steps  must  be  performed  to  quantita¬ 
tively  reconstruct  the  three-dimensional  transonic 
field  near  a  model  helicopter  blade  tip.  First, 
several  holographic  interferograms  must  be 
recorded  along  lines  perpendicular  to  the  rotor 
tip-path  plane  at  various  azimuthal  angles 
(fig.  3).  Data  must  then  be  extracted  from  the 
interferograms.  This  can  be  done  1)  by  hand, 

2)  by  using  a  graphic  tracing  tablet  [21],  or 

3)  by  using  a  system  that  digitizes  the  interfero¬ 
grams  and  extracts  fringe-order  numbers.  The 
digital  interferogram  evaluation  technique  was 
used,  and  is  presented  in  detail  by  Becker  and  Yu 
[22].  Finally,  the  data  are  transferred  as  input 
to  a  tomography  code,  which  computes  the  refrac¬ 
tive  index  at  specific  points  in  a  horizontal 
plane  above  the  blade  surface.  This  procedure  is 
repeated  in  several  planes  to  yield  a  reconstruc¬ 
tion  of  the  entire  three-dimensional  field. 

3. 1  Holographic  Interferogram  Recording 

The  holographic  system  for  recording  inter¬ 
ferograms  near  a  model  rotor  blade  was  assembled 
in  the  Aeroflightdynamios  Directorate  Anechoic 
Hover  Chamber.  Figure  4  shows  a  schematic  of 
the  optical  system  and  fig.  5  shows  the  Anechoic 
Hover  Chamber.  A  ruby  laser  with  a  20-nseo  pulse 
width,  a  694.3-nm  wavelength,  and  a  power  of  1  J 
"freezes"  the  rotating  blade  at  any  desired  azi¬ 
muthal  angle.  A  beam-splitter  divides  the  laser 
beam  into  two  separate  beams  at  the  laser  outlet. 

A  microscope  objective  lens  expands  the  object 
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beam  to  fill  a  61  cm-diam  spherical  mirror.  Since 
the  foci  of  both  the  objective  lens  and  the  spher¬ 
ical  mirror  coincide,  a  collimated  plane  wave 
forms  as  the  beam  passes  through  the  rotor  area. 
The  object  beam  then  strikes  a  second  61  cm-diam 
spherical  mirror,  emerges  as  a  converging  wave, 
and  illuminates  a  10.2  cm  by  12.7  cm  photographic 
plate.  The  reference  beam  is  lengthened  by 
causing  it  to  strike  several  plane  mirrors.  This 
beam  must  be  lengthened  so  that  the  difference  in 
the  path  lengths  of  the  object  and  reference  beams 
is  less  than  the  coherence  length  of  the  laser 
(one  of  the  very  few,  and  easily  met,  geometrical 
constraints  in  a  holographic  system).  The 
reference  beam  is  expanded  by  an  objective  lens, 
then  collimated  with  a  12.7  cm-diam  lens;  finally, 
it  is  directed  toward  the  film  so  that  it  overlaps 
the  object  beam. 

The  entire  procedure  is  conducted  from  out¬ 
side  the  hover  chamber,  once  the  optical  system  is 
aligned.  Firing  the  laser,  changing  the  photo¬ 
graphic  film  plates,  and  controlling  the  test 
conditions  are  all  done  by  remote  control.  Recall 
that  to  record  an  interferogram,  two  exposures  at 
different  times  (different  flow  states)  must  be 
made  on  a  single  film  plate.  The  film  records  the 
first  exposure  while  the  rotor  blade  remains  sta¬ 
tionary.  In  this  case,  the  air  has  no  velocity 
and  therefore  has  a  uniform  refractive-index  dis¬ 
tribution.  The  film  records  the  second  exposure 
while  the  blade  rotates  at  the  desired  speed.  The 
nonhomogeneous  refractive-index  distribution  in 
this  case  introduces  phase  changes  in  the  second 
object  wave,  producing  interference  fringes  on  the 
film  plate.  This  double-exposure  recording  pro¬ 
cedure  repeats  at  various  angles  around  the  flow 
by  synchronizing  the  laser  pulse  with  the  desired 
blade  position.  Because  of  the  long  optical  path- 
lengths  (27  m),  the  recording  system  is  very  sen¬ 
sitive  to  vibrations  of  the  optical  components. 

At  several  azimuthal  angles,  it  was  necessary  to 
record  multiple  interferograms  to  obtain  one  high- 
quality  interferogram.  The  photographic  plates 
are  then  removed  from  the  recording  system  in  the 
hover  chamber,  developed,  and  reconstructed  in  a 
laboratory  for  further  processing. 

Holographic  interferograms  record  the  flow 
near  a  hovering  1/7-soale  (geometric)  model  UH-'H 
rotor  with  untwisted  NACA  0012  airfoil  sections. 
The  blade  runs  at  a  tip  Mach  number  of  0.90  so 
that  the  flow  is  transonic  and  a  shock  wave  is 
present  [11.  This  model  normally  uses  two  blades; 
however,  in  views  along  the  span,  the  optical  beam 
would  pass  through  the  refract ive- index  field  of 
both  blades.  Because  tne  refractive-index  fields 
of  the  two  blades  are  inseparable  at  these  angles, 
a  single-bladed  ro' or  with  a  counterbalance  is 
used  instead  (fig.  6). 

Holographic  interferograms  near  a  tran¬ 
sonic  rotor  blade  are  sequentially  recorded  at 
40  different  viewing  angles.  Figure  7  shows  a 


simplified  schematic  plan  view  of  the  optical 
system  near  the  blade  area.  The  blade  rotates  in 
a  clockwise  direction  and  is  captured  at  any 
desired  viewing  angle  with  the  pulsed  laser.  The 
optical  system  translates  laterally  to  Insure  that 
the  laser  beam  passes  near  the  blade-tip  area. 

The  tomography  code  requires  flow  data  from 
certain  viewing  angles  within  a  180°  range. 
Numerical  simulation  results  [231  using  numerical 
computations  of  the  flow  [3]  suggest  recording 
interferograms  from  0=8°  to  0  =  40° ,  and 
from  0  =  140°  to  0  =  186°  in  2°  increments,  as 
defined  in  fig.  7.  The  missing  views,  0  =  42° 
to  0  I  138°,  were  presumed  to  have  very  few 
interference  fringes  and  were  not  utilized. 

Illustrated  in  fig.  7  are  examples  of  holo¬ 
graphic  interferograms  recorded  near  the  model 
blade  tip.  The  fringe  pattern's  appearance 
depends  on  the  viewing  direction.  Interferograms 
recorded  along  the  chord  (near  0  =  90°)  display 
very  few  interference  fringes,  since  the  optical 
rays  pass  through  the  field's  thinnest  (weakest) 
region.  No  observable  details  are  present  in 
these  views.  However,  in  views  along  the  span 
(near  0=0°  or  0  =  180°),  numerous  fringes  are 
visible,  because  the  optical  rays  pass  through  the 
longest  (strongest)  region  within  the  field.  The 
leading-edge  stagnation  point,  shock  structures, 
boundary-layer  separation,  and  wake  system  are 
clearly  visible.  In  particular,  a  lambda  shock 
(0  =  180°)  and  the  radiated  shock  (0  =  186°) 
appear  above  the  blade.  Several  interferograms 
are  described  in  detail  in  reference  24. 

3.2  Data  Extraction 

One  important  step  between  recording  the 
Interferograms  and  applying  tomographic  recon¬ 
struction  techniques  is  the  evaluation  of  the 
interferograms;  that  is,  reading  fringe  positions 
and  fringe  numbers.  Up  to  now,  most  interfero¬ 
grams  from  aerodynamic  tests  have  been  evaluated 
by  either  reading  fringe  numbers  and  their  posi¬ 
tions  manually  or  by  tracing  the  fringe  lines  by 
hand  with  the  help  of  a  tracking  device  (for 
instance  a  graphic  tablet).  However,  manual 
evaluation  is  a  very  time-consuming  and  inaccurate 
procedure.  Also,  it  is  evident  that  in  the 
current  application,  where  large  numbers  of  inter¬ 
ferograms  have  to  be  evaluated  at  several  hori¬ 
zontal  planes,  utilization  of  an  automatic  fringe¬ 
reading  procedure  would  enhance  the  evaluation  and 
would  make  the  interferometric  technique  a  much 
more  powerful  measurement  tool.  (See  reference  22 
for  a  detailed  discussion  of  this  procedure.) 

To  reconstruct  the  three-dimensional  flow 
field  from  these  interferograms,  using  tomographic 
reconstruction  techniques,  requires  the  fringe 
number-functions  along  cross  sections  of  a  plane 
pirallel  to  the  rotor  disk.  Data  from  all  the 
views  are  needed  to  reconstruct  the  refractive 


index  in  a  particular  plane  above  the  blade. 
Several  planes  have  to  be  reconstructed  to  get  a 
three-dimensional  flow-field  representation.  This 
procedure  requires  that  each  interferogram  be 
represented  in  a  form  that  allows  the  fringe-order 
function  to  be  computed  at  any  desired  point. 

The  evaluation  of  interferograms  using 
computer-aided  methods  can  be  subdivided  into 
the  following  steps:  >)  interferogram  digitizing 
and  image  enhancing,  2)  fringe  segmenting  and 
fringe-coordinate  extracting,  3)  obtaining  fringe 
fields  from  several  magnified  views  and  merging, 

4)  fringe  numbering  and  correcting  fringe  discon¬ 
nections,  5)  transforming  coordinates,  6)  inter¬ 
polating  and  extrapolating  fringe-number 
functions.  In  the  following,  a  system  Is  briefly 
described  for  evaluating  interferograms  that 
incorporates  the  evaluation  steps  mentioned  above. 

3.2.1  Hardware  Components.  An  image-processing 
system  (De  Anza  IP-6400)  connected  to  a  VAX  11/780 
host  computer  provides  the  main  hardware  necessary 
for  digitizing  the  interferograms  and  for  doing 
some  image-enhancement  processing.  The  resolution 
of  the  system  is  512  .  512  pixels  with  an  8-bit 
intensity  range.  Currently  it  is  equipped  with 
two  of  possible  four  memory  planes,  as  well  as 
with  a  graphic  and  an  alphanumeric  overlay.  A 
frame-grabbing  unit  can  digitize  a  frame  of  a 
Video  signal  in  real  time  (1/30  sec).  A  black  and 
white  video  camera  (MTl  series  68)  with  a 
resolution  (bandwidth)  of  18  MHz  is  connected  to 
this  Input  channel . 

The  system  has  an  arithmetic-logic  unit 
(ALU),  with  which  real-time  addition,  subtraction, 
or  comparison  of  one  or  more  image  planes  may  be 
made.  The  contents  of  each  memory  plane  may  be 
routed  through  lookup  tables  before  being  input  to 
the  ALU,  to  the  video-output  processor,  or  to 
another  plane.  The  actual  contents  of  any  image 
plane,  or  of  a  combination  of  image  planes,  is 
output  via  a  video-output  processor  and  can  be 
shown  on  a  color-display.  Each  channel  has  its 
own  color-mapping  tables.  A  Joystick  control 
device  is  used  for  interactive  input.  It  controls 
two  cursors,  which  may  be  used  in  a  number  of 
different  operating  modes.  A  color  print  system 
(Dunn  Instruments  model  631)  serves  as  a  hardcopy 
device  for  the  color  monitor. 

3.2.2  Digitization  and  Processing.  During  the 
recording  of  the  holographic  interferograms,  two 
fiducial  points  were  marked  in  the  image  plane 
with  a  known  position  relative  to  the  rotor 
system.  This  allowed  a  subsequent  coordinate 
tr.ansformat ion  i-to  a  space-fixed  coordinate 
system.  E.ich  interferogram  was  then  digitized, 
and  several  subsections  were  enlarged,  depending 
or.  the  fringe  density  in  the  interferograms.  A 
sc.iie.  aligned  to  the  fiducial  marks,  was  always 
digit. zed  together  with  the  interferogram  so  that 


the  position  and  magnificaticn  of  the  enlarged 
segments  could  be  identified.  Before  segmentation 
and  fringe-coordinate  extraction  are  applied,  it 
may  be  useful  to  do  some  image  enhancement  to 
reduce  noise.  Frame-averaging  methods  increase 
the  signal-to-noise  ratio  of  the  imaging  elec¬ 
tronic  components,  thereby  Improving  the  picture 
quality. 

3.2.3  One-Dimensional  Evaluation.  One- 
dimensional  interactive  processing  may  be 
employed,  if  only  one  two-dimens io.nal  plane  of  the 
flow  field  is  reconstructed;  hence,  only  one  cross 
section  through  each  interferogram  has  to  be 
known.  In  this  case  some  of  the  expense  required 
for  the  full  two-dimensional  evaluation  of  the 
fringe  pattern  may  be  avoided.  But  the  procedure 
has  to  be  interactive,  because  Iccal  information 
(l.e.,  knowledge  of  the  fringe  locations  along  one 
line)  is  not  sufficient  to  number  the  fringes 
correctly  or  to  detect  lost  fringes  or  other  dis¬ 
tortions.  In  an  interactive  procedure,  informa¬ 
tion  such  as  acceptance  of  the  segmentation  and 
assignment  of  fringe  numbers  has  to  be  given  by 
the  user. 

A  program  featuring  the  one-dimensional  eval¬ 
uation  digitizes  and  preprocesses  an  interferogram 
and  does  a  fringe  segmentation  along  a  line,  or 
along  a  set  of  lines,  through  the  field  of  view. 
The  result  of  the  segmentation  procedure  (a  binary 
fringe  pattern)  is  written  pack  to  the  image 
screen  for  monitoring  reasons  (see  figs.  8a  and 
8b).  Users  may  interactively  change  the  accep¬ 
tance  threshold  or  edit  the  segmented  cross  sec¬ 
tion  in  heavily  distorted  regions  before  they 
continue  to  the  numbering  section.  The  cursor  can 
be  moved  to  each  part  of  the  segmented  line,  and 
fringe  numbers  may  be  assigned  interactively  by 
various  commands.  The  task  assigning  fringe  num¬ 
bers  is  supported  by  color  coding  the  black  parts 
of  the  fringes;  this  shows  whether  the  fringe- 
order  function  increases  or  decreases  by  1,  or  if 
a  discontinuity  is  present  between  adjacent 
fringes.  The  resulting  fringe-order  function  may 
also  be  plotted  onto  a  graphic  terminal.  A  post¬ 
processing  program  merges  all  the  data  taken  from 
different  enlarged  portions  of  an  interferogram 
(as  shown  in  figs.  8a  and  8b)  and  uses  a  spline- 
approximation  to  interpolate  between  the  fringes. 
The  final  result  is  extracting  the  fringe  order 
function  at  a  desired  height  above  the  blade  sur¬ 
face  as  shown  in  figure  8c. 


3. 2. A  Two-Dimens ional  Interactive  Processing. 
Often  the  knowledge  of  only  a  cross  section  across 
the  field  of  view  is  not  sufficient  for  inter¬ 
preting  a  given  flow  problem;  instead,  the  fringe 
patterns  have  to  be  evaluated  over  the  whole  field 
of  view.  To  facilitate  the  two-dimensional  inter¬ 
ferogram  evaluation,  additional  algorithms  have 
been  implemented  to  segment  fringes,  to  extract 


each  fringe's  coordinates  (trace  the  fringe 
sides),  and  to  represent  the  whole  fringe  field  as 
a  polygonal  data  structure.  Methods  for  numbering 
these  line  fields  and  detecting  extracting  errors 
have  been  developed.  With  a  two-dimensional 
Interpolation  scheme,  fractional  fringe-order 
numbers  may  then  be  estimated  at  any  given  point 
in  the  field  of  view. 

Figure  9  shows  an  example  of  a  typical 
fringe-coordinate  extraction.  An  interferogram 
(fig.  9a)  is  digitized  over  the  whole  field  of 
view  with  low  resolution.  The  white  lines  over¬ 
laying  the  display  represent  the  extracted  fringe 
sides,  and  the  black  line  represents  the  border 
polygon.  The  border  polygon  was  input  by  hand 
before  the  coordinate  extraction  in  order  to  mask 
off  background  objects  and  the  part  of  the  fringe 
field  with  high  fringe  density.  (In  our  implemen¬ 
tation,  all  the  graphics  are  in  color  and  there¬ 
fore  easier  to  distinguish  from  the  displayed 
image. ) 

As  was  discussed  before,  it  may  be  necessary 
to  digitize  enlarged  subsections  to  resolve 
regions  with  high  fringe  densities.  The  resulting 
fringe  polygonal  field  la  then  combined  out  of  all 
the  individual  digitized  and  processed  sub¬ 
sections,  starting  with  the  one  with  highest  reso¬ 
lution.  To  handle  the  problem  of  partly  overlap¬ 
ping  polygonal  line  fields,  a  boundary  polygon  is 
maintained  for  each  subsection.  This  defines  the 
definition  area  for  each  fringe  polygonal  field. 
Upon  combination  of  two  overlapping  fields,  the 
border  polygon  of  the  one  with  higher  resolution 
(priority)  is  used  to  intersect  the  lines  of  the 
other  field.  Both  line  fields  are  then  connected 
at  those  intersection  points.  The  new  border 
polygon  is  the  border  of  the  combined  areas. 

Figure  9b  shows  an  enlarged  digitized  part  of  the 
interferogram  shown  in  fig.  9a,  with  the  extracted 
fringe  polygons  overlaying. 

3.3  CAT  Reconstruction 

Fringe-order  functions  are  transferred  to  the 
filtered  back-projection  CAT  code  which  computes 
the  refractive-index  field  at  specific  points  in  a 
chosen  horizontal  plane  above  the  blade  surface. 
The  code  assumes  refractionless  light  rays;  there¬ 
fore,  each  horizontal  plane  can  be  treated  inde¬ 
pendently,  even  though  data  for  each  plane  is 
taken  from  one  set  of  interferograms .  The  pres¬ 
sure  coefficient  is  computed  by  first  converting 
refractive  index  to  density,  using  the  Giadston- 
Dale  relation: 

o/pQ  :  (n  -  '  )/(n|j  -  1 ) 

Density  is  then  converted  to  pressure  coefficient 
from  a  form  of  Bernoulli's  equation  for  steady 
(with  respect  to  the  rotation  blade),  compres¬ 
sible,  isentropic  flow.  The  procedure  is  repeated 


in  several  planes  above  the  blade  to  reconstruct 
the  entire  three-dimensional  field  near  the  axjdel 
blade  tip. 

4.  RECONSTRUCTION  RESULTS 

The  holography-CAT  reconstruction  of  the 
blade-tip  pressure  coefficient  (Cp)  field 
is  compared  with  numerical  computations  and  laser 
velocimeter  (LDV)  measurements.  The  computa¬ 
tions  used  here  are  conservative,  mixed-difference 
solutions  of  the  transonic  small-disturbance 
equation  [31.  The  LDV  measurements  obtained  by 
Owens  et  al.  [25)  were  recorded  using  the  same 
test  facility  and  blades  as  the  holographic 
experiment.  All  results  were  obtained  with  a  tip 
Mach  number  of  0.90.  Results  from  these  three 
sources  are  presented  in  three  horizontal  planes 
above  the  airfoil  center  axis.  Three  plot  types 
are  used  to  visualize  the  flow  field.  First,  Cp 
contours  are  displayed  comparing  the  holography- 
CAT  reconstruction  and  the  numerical  computa¬ 
tions.  The  contours  are  given  in  plan  view,  where 
the  blade's  leading  and  trailing  edges  are 
X/C  =  0.0  and  X/C  =  1.0,  respectively.  The 
blade  tip  is  located  at  =  1.0,  the  rotation 

center  is  at  R/R^  =  0.0,  and  the  blade  rotates  in 
a  clockwise  direction.  Second,  perspective  views 
are  displayed  in  which  negative  Cp  values  are 
plotted  along  the  vertical  axis.  The  data  and 
geometry  are  identical  to  the  contour  plots, 
though  the  data  are  viewed  from  near  the  rotor 
hub.  Third,  Cp  distributions  from  all  three 
sources  are  shown  at  several  radial  locations. 

Figure  10  compares  Cp  contours  derived  from 
the  holographic-CAT  method  and  the  numerical  com¬ 
putations  near  the  blade  surface  (Y/C  =  0.08). 

Both  results  display  high-Cp  regions  near  the 
leading  and  trailing  edges.  They  also  display  a 
low-Cp  region  over  the  blade  surface,  containing  a 
shock  at  approximately  X/C  i  0.60  near  the  blade 
tip.  The  general  contour  shapes  show  strong 
resemblance  except  near  the  blade  tip  (roughly  the 
last  5t  of  blade  span).  The  minimum  Cp  region 
appears  at  the  shock  foot  in  the  computational 
analysis,  but  it  appears  closer  to  the  leading 
edge  and  farther  from  the  shock  foot  in  the 
holography-CAT  results. 

Figure  11  shows  C  values  for  the  same  plane 
in  perspective  view.  Again,  the  general  flow 
shapes  appear  very  similar.  The  major  difference 
here  between  the  two  results  is  the  roughness 
(minor  "ridges")  in  the  reconstructed  flow.  This 
may  be  due  to  reconstruction  artifacts  caused  by 
noncontinuous  data  (interferograms  recorded  in 
2  azimuthal  increments)  or,  more  likely,  by  noise 
(erroneous  fringes)  in  the  interferogram  data 
caused  by  optical  component  motion.  Figure  12 
compares  C  distributions  from  the  holographic-CAT 
reconstruction,  numerical  computations,  and  LDV 
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measurements  (where  available).  The  roughness  of 
the  reconstructed  flow  can  be  seen  throughout  the 
figures.  Also,  the  holography -CAT  method  deter¬ 
mines  the  shock  location  to  be  slightly  aft  (3-5%) 
of  the  numerical  results  or  LDV  measurements  (see 
fig.  12c).  The  major  differences  can  be  observed 
in  fig.  I2d,  where  the  discrepancies  at  the  lead¬ 
ing  edge  and  over  the  blade  surface  are  clearly 
visible.  Note  that  all  three  methods  give  dif¬ 
ferent  results,  especially  over  the  blade  surface. 

Figure  13  compares  Cp  contours  from  both  the 
holography-CAT  method  and  the  numerical 
computations  at  Y/C  :  0.22  above  the  blade 
center  axis.  The  general  Cp  contours  show  an 
excellent  agreement  in  both  shape  and  magnitude 
throughout  this  plane.  The  minimum  C  region 
above  the  blade  (near  X/C  =  0.50,  R/Bq  =  0.96) 
matches  much  closer  in  this  plane  than  in  the 
plane  near  the  blade  surface  (fig.  10).  The 
perspective  view  of  fig.  19  also  shows  an 
excellent  agreement  in  shape  and  magnitude. 

Again,  the  most  noticeable  difference  in  the 
reconstructed  flow  is  the  extra  ridges  at  the  same 
locations  and  orientations  as  seen  in  the  previous 
plane  (fig.  11a).  The  holography-CAT  Cp 
distributions  in  this  plane  (fig.  15)  compare 
favorably,  especially  inboard  the  blade  tip  (see 
figs.  15a-15c).  Even  at  the  blade  tip  (fig.  15d), 
the  holography-CAT  reconstruction  shape  and  magni¬ 
tude  agree  more  closely  with  the  other  methods 
than  in  fig.  12d. 

Finally,  fig.  16  compares  Cp  contours  from 
the  holography-CAT  method  and  numerical  computa¬ 
tions  for  Y/C  :  0.99  above  the  blade.  The  Cp 
contour  shapes  and  magnitudes  are  similar  in  all 
regions,  except  that  the  holography-CAT  results 
show  the  minimum  Cp  point  to  be  slightly  (1%  span) 
outboard  of  the  numerical  result.  The  extra 
ridges  seen  in  perspective  view  (fig.  17a)  are  in 
approximately  the  same  location  and  orientation  as 
in  previous  planes;  however,  the  magnitude  of  the 
ridges  has  decreased.  Figure  18  shows  Cp  distri¬ 
butions  for  Y/C  =  0.99.  There  is  good  agreement 
at  all  locations.  The  most  noticeable  differences 
are  the  shock  location  in  fig.  I8a  (holography -CAT 
results  slightly  aft)  and  the  magnitude  in 
fig.  l8d  (holography-CAT  results  slightly  larger). 

Overall,  the  comparison  of  the  holography-CAT 
reconstruction  with  numerical  computations  and  LDV 
measurements  is  extremely  encouraging  in  most 
regions  of  the  flow  field.  However,  comparisons 
with  other  numerical  code  computations  and  other 
experimental  data  (i.e.,  hot  wire  measurements) 
are  required  before  a  final  evaluation  can  be  made 
about  the  holography-CAT  results.  Several  dis¬ 
crepancies  must  be  resolved.  First,  the  shock 
location  must  be  confirmed  (perhaps  through 
pressure- instrumented  blade  measurements). 

Second,  the  extra  ridges  appearing  in  the  recon¬ 
structed  results  must  be  reduced  or  eliminated. 

The  ridges  are  most  likely  caused  by  poor  data 


recorded  in  the  interferograms.  To  improve  the 
interferograms'  quality,  a  modified  optical  system 
(for  both  hover  and  forward  flight  testing)  is 
necessary.  In  addition,  all  methods  were  unable 
to  predict  the  expected  shock  strength  beyond  the 
blade  tip.  Acoustic  measurements  indicate  a  much 
stronger  radiated  shock  than  is  indicated  by  these 
three  results. 

5.  CONCLUDING  REMARKS 

The  holographic  Interferometry  computer- 
assisted  tomography  technique  proved  to  be  a 
highly-effective  way  of  measuring  the  three- 
dimensional,  transonic  flow  field  near  a  model 
rotor-blade  tip.  Results  from  this  method  compare 
favorably  with  those  of  numerical  computations  and 
laser  velocimeter  measurement  except  very  near  the 
tip  region.  In  other  regions,  the  Cp  distribu¬ 
tions  along  the  chord  are  similar  in  both  shape 
and  magnitude.  However,  the  results  from  the 
technique  must  be  further  verified  against  other 
experimental  data. 

Since  this  is  the  first  successful  imple¬ 
mentation  of  the  holographic  interferometry, 
computer-assisted  tomography  method  in  rotor  flow 
studies,  many  improvements  are  indicated.  For 
example,  the  optical  system  must  be  improved  so 
that  better  quality  interferograms  can  be 
recorded,  and  an  automatic  fringe-reading  tech¬ 
nique  must  be  completed  so  that  the  time  required 
to  evaluate  interferograms  can  be  shortened.  Upon 
verification  of  these  results  and  after  the  system 
is  improved,  measurements  of  other  model  rotor- 
blade  flow  fields,  including  those  of  forward 
flight,  will  be  performed. 
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(a)  (b) 

Fig.  1  Back-projection,  a)  One  projection  of  an  absorbing  disk;  b)  back-projecting 
consists  of  smearing  each  projection  back  along  the  direction  in  which  the  original 
projection  was  made  (from  ref.  19). 
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Fig.  2  Filtered  back  projection.  a)  The  projection  data  are  convolved  (filtered) 
with  a  suitable  processing  function  before  back-projection;  b)  three  back- 
projections  of  an  absorbing  disk  (from  ref.  19). 


Fig.  11  Pressure-coefficient  values  for  Y/C  :  0.08  above  blade  centerline. 
a)  Holography-CAT  reconstruction;  b)  numerical  computation. 
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Fig.  12  Pressure-coefficient  distributions  at  six  radial  locations  for 

Y/C  =  0.08  above  blade  centerline,  a)  R/R^  -  0.90;  b)  R/Rq  =  0.92;  c)  R/R^  =  0.95; 

d)  R/Rq  =  1.00;  e)  R/R.  =  1.05;  f)  R/R^  =  1.10. 


Fig.  13  Pressure-coeficient  contours  in  plan  view  for  Y/C  =  0.22  above  blade 
centerline:  contour  interval  -  0.1.  a)  Holography-CAT  reconstruction;  b)  numerical 
computation. 


Fig.  14  Pressure-coefficient  values  for  Y/C  =  0.22  above  blade  centerline 
a)  Holography-CAT  reconstruction;  b)  numerical  computation. 
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Fig.  16  Pressure-coefficient  contours  in  plan  view  for  Y/C  =  0.i(9  above  blade 
centerline:  contour  interval  *  0.07.  a)  Holography-CAT  reconstruction; 
b)  numerical  computation. 


Fig.  17  Pressure-coefficient  values  for  Y/C  =  0.49  above  blade  centerline, 
a)  Holography-CAT  reconstruction;  b)  numerical  computation. 
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Fig.  18  Pressure-coefficient  distributions  at  six  radial  locitions  for 
Y'C  =  0.49  above  blade  centerline,  a)  =  0.90;  b)  R  =  0.9.?:  c)  R  R^  = 

d)  R/R  =  1.00;  e)  R/R  =  1.05:  f)  R- R„  =  1.10. 
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